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LONG-TERM  GOALS 

The  long-term  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world's  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water.  An 
associated  goal  is  the  understanding  of  the  absorption  and  backscattering  properties  of  marine  particles 
in  terms  of  the  distributions  of  their  size,  shape,  and  composition. 

OBJECTIVES 

The  inherent  optical  properties  (lOPs)  of  marine  particles  are  most-often  modeled  as  homogeneous 
spheres  using  Mie  Theory.  Although  this  approach  has  been  fruitful,  the  next  logical  step  in  modeling 
marine  particles  is  to  abandon  the  normally-employed  spherical  approximation  and  use  more  realistic 
approximations  to  their  shape.  The  advent  of  computer  codes  capable  of  handling  more  complex 
shapes,  and  the  increased  computational  speeds  now  available,  suggest  that  particle  modeling 
employing  simple  non-spherical  shapes,  e.g.,  disks,  rods,  etc.,  could  become  routine.  For  example, 
Gordon  and  Du  (2001)  used  a  two-disk  model  to  try  to  reproduce  the  backscattering  by  coccoliths 
detached  from  E.  huxleyi;  however,  they  found  that,  while  the  resulting  spectral  variation  of  the 
backscattering  cross  section  agreed  with  experiment,  its  magnitude  was  low  by  a  factor  of  2-3.  Such 
simple  shapes  are  still  at  best  poor  approximations  to  real  particles,  as  real  particles  display  significant 
inhomogeneity. 

It  is  often  assumed  that  inhomogeneous  particles  can  be  replaced  by  homogeneous  particles  along  with 
employing  an  “effective”  refractive  index;  however,  my  study  of  thin  disks  showed  the  importance  of 
periodic  angular  fine  scale  structure  and  demonstrated  that  if  the  inhomogeneities  are  too  large,  the 
effective  refractive  index  assumption  is  invalid  and  the  backscattering  increases  significantly  over  that 
expected  based  on  the  effective  index  hypothesis.  Thus,  in  attempting  to  model  with  shapes  other  than 
spherical,  we  are  still  left  with  the  question  of  the  relative  importance  of  gross  morphology  and  fine 
structure.  To  apply,  with  confidence,  scattering  computations  based  on  simple  particle  shapes  to 
explain  or  interpret  observed  particle  backscattering,  it  is  important  to  understand  how  deviations  from 
simple  shapes  or  how  non-homogeneity  of  composition,  affect  the  scattering  and  absorption  of  light. 
The  objective  of  this  study  is  to  continue  the  development  of  such  an  understanding. 
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APPROACH 


I  use  detached  coccoliths  from  the  coccolithophorid  huxleyi  (Figure  1)  as  a  case  study  for  applying 
non-spherical  shapes  to  the  computation  of  backscattering.  This  particle  was  chosen  because  (1)  its 
composition  is  known  and  homogeneous  (refractive  index  m  «  1 .20),  (2)  its  gross  shape  is  known  (it 
resembles  a  disk  or  two  parallel  disks),  (3)  its  optical  properties  are  known,  and  (4)  it  has  a  complex, 
quasi-periodic  internal  structure  on  a  sub-visible-wavelength  scale.  Gordon  (2006)  focused  on  the 
quasi-periodic  structure  of  the  top  disk  (in  Figure  1)  and  on  the  influence  of  the  curvature  of  the  lower 
plate.  (The  influence  of  the  curvature  was  found  to  be  small  and  will  not  be  discussed  here.) 


Figure  1:  SEM  image  of  coccolith  a  detached  from  E.  huxleyi  showing 
a  structure  grossly  resembling  two  parallel  disk-like  plates. 

The  upper  plate  resembles  a  wheel  with  spokes. 

The  spoke  separation  is  50  to  150  nm. 


I  use  the  term  “gross  morphology”  to  indicate  a  smooth  homogeneous  particle  having  approximately 
the  same  overall  shape  as  the  particle  in  question  (e.g.,  a  single  disk  or  two  parallel  disks  as  a  model 
for  a  detached  coccolith).  I  use  “fine-scale  structure”  to  indicate  deviations  from  the  gross  morphology 
(e.g.,  the  coccolith's  periodic  radial  structures  resembling  the  openings  between  the  spokes  of  a  wagon 
wheel).  My  goal  is  to  understand  how  the  fine-scale  structure  can  induce  deviations  in  the 
backscattering  characteristic  of  a  given  gross  morphology.  Since  the  physical  example  of  interest  is  a 
coccolith,  I  limit  the  examination  to  particles  with  the  gross  morphology  of  a  disk. 

To  investigate  the  influence  of  periodic  fine  structure  in  a  disk-like  object  on  backscattering,  I  started 
with  a  homogeneous  disk  and  removed  sectors.  Specifically,  the  disk  was  divided  into  equal  angle 
sectors  of  angle  Aa,  and  alternate  sectors  were  removed.  The  angle  A«  was  given  by  A«  =2iitT , 
where  n  is  an  integer.  I  refer  to  these  objects  as  “pinwheels.”  (See  Figure  2,  ^  =  0  column,  for  an 
example  of  a  pinwheel  with  n  =  5.)  If  we  let  s  be  the  arc  length  of  the  open  (or  closed)  regions  at  the 
perimeter  of  the  pinwheel,  then  s  =  DAa  12,  where  D  is  the  diameter 
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^  =  0 


£  =  0.5 


1.0 


Figure  2:  The  individual  rows  provide  four  realizations  of  the  aperiodic  pinwheels 
for  n  =  5  and  various  values  of  e.  Periodic  pinwheels  are  shown  in  the  first  column. 
For  the  aperiodic  pinwheels,  the  individual  vanes  have  random  angular  extent. 
The  departure  from  the  periodic  pinwheel  increases  with  increasing  & 


of  the  disk  (1.5  |u,m).  The  values  of  s  for  the  various  oases  that  1  examined  were  suoh  that  at  a 
wavelength  (A)  of  400  nm  in  vaouum  (300  nm  in  water),  as  n  progresses  from  4  to  7,  5  takes  on  the 
values  A,  A/2,  A/4,  and  A/%  in  water.  One  of  the  main  goals  of  my  study  was  to  determine  if  a 
relationship  exists  between  s  and  A  where  the  periodio  struoture  beoomes  important  (or  unimportant)  to 
the  baoksoattering.  1  oomputed  the  baoksoattering  oross  seotion  o*  of  pinwheels  using  the  disorete- 
dipole  approximation  (Draine,1988;  Draine  and  Flatau,  1994).  The  oomputations  showed  that  their 
baoksoattering  oross  seotion  was  nearly  identioal  to  that  of  a  homogeneous  disk  of  similar  size  (but 
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with  m  reduced  from  1.20  to  meff=  1.10)  as  long  as  V/lwater  <  0.25.  Pinwheels  satisfying  this  criterion 
also  scattered  in  the  same  manner  as  disks  having  half  of  the  mass  removed  from  random  locations 
within  the  disk.  Thus,  when  V/lwater  <0.25,  there  is  essentially  no  difference  in  backscattering  between 
the  periodic  structure  and  a  structure  of  small  random  voids  (in  disk-like  particles)  with  the  same  total 
mass.  In  this  regime  the  backscattering  is  totally  governed  by  the  particle’s  gross  morphology  and 
effective  index  (meff)-  For  V/lwater  >  0.25  departures  from  a  homogeneous  disk  are  observed  and 
manifest  as  a  significant  increase  (many  times)  in  backscattering:  the  fine-scale  structure  is  as 
important  as  the  gross  morphology  in  this  regime. 

In  the  present  work  1  continue  exploration  of  the  backscattering  by  objects  possessing  small-scale 
structure.  1  focus  on  (1)  extending  the  sectorized  disk  studies  (Gordon  2006)  to  larger  t/Zwater  in  order 
to  understand  how  large  cr*  will  eventually  become  when  V/lwater  >0.25,  (2)  extending  the 
computations  to  sectorized  disks  with  aperiodic  structure  to  try  to  understand  the  influence  of  non- 
uniform  spacing  of  the  spokes  in  the  upper  plate  of  the  coccolith  (Figure  1)  and  (3)  building  a  fine- 
scale  model  of  a  detached  coccolith  for  computation  of  the  coccolith-specific  cr*. 

WORK  COMPLETED 

1  formed  aperiodic  pinwheels  through  a  perturbation  of  the  purely  periodic  pinwheel  effected  in  the 
following  manner.  First,  as  before,  the  disk  is  divided  into  purely  periodic  sectors,  the  angular 
boundaries  of  which  are  designated  by  the  2”  angles  ap.  The  individual  boundary  angles  are  then 
perturbed  to  ai  according  to 
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where  0  <  s<  1  is  a  constant  and  -1/2  <  p<  1/2  is  a  random  number  with  a  uniform  probability 
density.  Then,  the  material  of  the  disk  is  removed  from  ever  other  sector,  yielding  a  pinwheel  with  a 
quasi-periodic  structure.  Four  realizations  (each  based  on  a  difference  sequence  of  pseudorandom 
numbers)  of  such  pinwheels  for  n  =  5  are  provided  in  Figure  2  for  0.5  and  1.0  in  the  second  and 
third  columns,  respectively.  Defining  07  to  be  the  standard  deviation  in  the  angle  a/,  we  find  07  =  12 
Aa  «  0.3f  Aa,  where  Aa  =  lull".  Likewise  defining  (lAa  to  be  the  standard  deviation  of  the  removed 
(or  occupied)  sector  angles,  <JAa  =  t:^a  «  0.4^  Aa.  Thus,  for  e  =  0.5  and  1.0,  the  relative  standard 
deviation  in  angle  of  the  removed  (or  occupied)  sectors  is  20  and  40%,  respectively.  1  examined  three 
aperiodic  pinwheels.  The  first  has  a  diameter  (D)  of  1.50  pm,  a  thickness  (t)  of  0.15  pm,  and  n  =  5 
(Figure  2),  the  second  has  D  =  2.75  pm,  t  =  0.05  pm,  and  n  =  6,  and  the  third  has  £>=1.50  pm,  t  =  0.05 
pm,  and  n  =  5.  The  larger-diameter  disk  is  similar  in  size  to  the  distal  shield  (top  plate)  of  individual 
E.  huxleyi  coccoliths  (Figure  1). 

The  individual  realizations  are  labeled  by  the  one  minus  the  position  in  a  string  of  pseudorandom 
numbers  where  the  sampling  for  begins.  Thus,  for  realization  0000  the  sampling  begins  with  the 
first  number,  for  realization  1000  it  begins  with  the  1001**'  number,  etc.  This  method  of  creating 
aperiodic  pinwheels  does  not  yield  structures  with  the  same  volume  (mass)  as  the  associated  periodic 
pinwheel.  The  variation  of  the  volume  for  a  given  ^  can  be  as  much  as  25%  for  the  smaller  (Z)  =  1.5 
pm,  n  =  5)  and  7%  for  the  larger  {D  =  2.75  pm,  n  =  6)  pinwheels.  The  reduction  in  dispersion  from 
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the  smaller  to  the  larger  is  due  to  the  inerease  in  n,  whieh  doubles  the  number  of  seetors,  inereasing  the 
probability  that  the  individual  realizations  have  a  volume  eloser  to  the  mean. 

The  scattering  computations  were  carried  out  using  the  discrete-dipole  approximation  (DDA)  operating 
on  a  40  CPU  cluster.  The  accuracy  of  the  DDA  for  randomly  oriented  particles  is  governed  by  two 
issues:  (1)  employing  a  sufficient  number  of  dipoles  to  solve  the  electromagnetic  scattering  problem 
for  a  given  orientation;  and  (2)  employing  a  sufficient  number  of  orientations  for  performing  the 
orientational  average.  A  measure  of  the  number  of  dipoles  is  related  to  d,  the  spacing  between  the 
dipoles.  One  wants  d  to  be  substantially  smaller  than  the  wavelength.  The  smaller  d,  the  more  dipoles 
are  required  to  fill  the  volume  of  the  particle.  A  convenient  measure  of  the  spacing  in  regard  to  the 
wavelength  is  \m\kd,  where  m  is  the  refractive  index  and  k  =  In! A.  Gordon  and  Du  (2001)  showed 
that,  for  a  homogeneous  disk  with  D  =  2.1  pm,  using  -'5000  orientations  for  the  orientational  average, 
the  error  in  the  backscattering  cross  section  (o*)  was  of  the  order  of  5%  for  \m\kd  <  0.5,  and  decreased 
rapidly  for  smaller  values  of  \m\kd.  In  the  present  work,  I  always  used  enough  dipoles  to  keep  \m\kd  < 
0.5  (and  often  <  0.4)  and  a  number  of  orientations  that  would  provide  approximately  the  same 
averaging  accuracy  as  the  uniform  disk  with  D  =  2.7  pm.  For  the  aperiodic  pinwheels  extensive  testing 
showed  that  90,000  orientations  were  sufficient  to  perform  the  averaging. 

RESULTS 

Although  I  carried  out  the  computations  for  pinwheels  of  three  sizes,  I  will  only  discuss  the  first  {D  = 
1.50  pm,  t=  0.15  pm,  and  n  =  5)  in  detail  here.  Figure  3  provides  the  results  of  the  backscattering 
cross  section  computations.  Also  shown  for  comparison  is  the  result  for  a  homogeneous  disk  of  the 
same  size  {D  and  t).  Note  that  the  homogeneous  disk  has  twice  the  volume  (mass)  of  the  periodic 
pinwheel  {s=  0)  and  approximately  twice  the  volume  of  the  aperiodic  pinwheels  {s  >  0). 

Consider  first  the  periodic  pinwheels.  For  these,  5  =  XI A  occurs  when  t/X  =  0.255.  For  t/X  larger  than 
this  value,  the  backscattering  increases  rapidly  with  decreasing  X  and  then  undergoes  a  series  of 
maxima  and  minima  with  progressively  increasing  backscattering  at  each  maximum.  The 
backscattering  at  the  maxima  is,  in  magnitude,  approximately  that  at  the  maxima  for  the  homogeneous 
disk  (twice  the  volume  or  mass  of  the  pinwheel).  These  maxima  are  the  result  of  interference  of  the 
fields  scattered  by  the  individual  vanes  of  the  pinwheel  as  they  occur  in  the  Rayleigh-Gans 
approximation  as  well  (although  only  the  first  maximum  occurs  at  the  same  position  [Gordon,  2007]). 
In  the  case  of  the  smaller  aperiodic  pinwheel,  when  s=  0.5,  its  backscattering  closely  follows  the 
periodic  pinwheel,  with  the  dispersion  of  backscattering  reaching  20%  at  the  smallest  wavelength, 
while  when  ^  =  1.0  the  dispersion  is  somewhat  larger.  It  is  interesting  to  note  that  near  the  position  of 
the  first  maximum,  the  aperiodic  pinwheel  backscatters  less  than  the  periodic  pinwheel,  and  the 
deviation  between  the  two  increases  with  increasing  s.  This  behavior  would  be  expected  under  the 
hypothesis  that  the  maxima  in  the  periodic  case  results  from  constructive  interference  of  light 
interacting  with  the  individual  vanes  of  the  pinwheel  -  when  the  spacing  and  angular  size  of  the  vanes 
becomes  random  the  constructive  interference  is  reduced. 
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Figure  3:  The  backscattering  cross  section  of four  realizations  of  the  aperiodic  pinwheels  shown  in 
Figure  1  (D  =  1.5  pm)  compared  with  that  for  a  periodic  pinwheel  and  for  a  homogeneous  disk  of 
the  same  size  as  a  function  of  the  thickness  divided  by  the  wavelength.  The  left  panel  is  for  £=0.5 
and  the  right  panel  is  for  £  =  1.0.  The  periodic  pinwheel  shows  increasing  (and  oscillatory) 
backscattering  with  decreasing  wavelength.  The  individual  realizations  of  aperiodic  pinwheels  all 
roughly  follow  their  periodic  counterpart,  but  with  dispersion  that  increases  with  increasing  & 


Figure  4:  The  backscattering  cross  section  {at)  of  an  equal-number  mixture  of  the  four 
realizations  of  the  aperiodic  pinwheels  compared  with  that  for  a  periodic  pinwheel  (£=0).  This 
shows  that  the  £=0.5  cases  are  close  to  the  periodic  pinwheel,  but  for  £=1.0,  the  deviation  is 
considerable,  especially  near  the  first  maximum  for  £=0. 
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It  is  important  to  understand  that  real  biologieal  partieles  (e.g.,  E  huxleyi  eoecoliths)  would  display 
similar  dispersions  in  baekseattering.  In  faet,  if  pinwheels  were  to  represent  real  biological  particles, 
samples  would  be  expected  to  consist  of  a  number  of  realizations  of  their  aperiodicity.  In  this  regard, 
the  average  cr*  (denoted  by  (o*))  is  more  important  than  that  for  any  given  realization.  Figure  4 
compares  the  (o*)  for  the  four  realizations  of  the  aperiodicity  examined  here  with  the  associated 
periodic  pinwheel.  It  clearly  shows  that  the  main  difference  between  (cr*)  for  the  small  periodic  and 
aperiodic  pinwheels  (left  panel)  occurs  near  the  maxima  in  the  baekseattering,  and  that  near  the  first 
maximum  (but  not  the  second)  the  difference  increases  as  deviation  from  periodicity  increases  (i.e.,  as 
f  increases). 

The  computations  for  the  two  other  examples  of  aperiodic  pinwheels  {D  =  1.5  pm,  t  =  0.05  pm,  n  =  5, 
and  D  =  2.75  pm,  t  =  0.05  pm,  n  =  6)  revealed  that  as  DIt  becomes  larger,  the  dispersion  among 
realizations  decreases  somewhat,  and  (cr*)  becomes  very  close  to  that  of  the  associated  periodic 
pinwheel. 

IMPACT/APPLICATIONS 

I  examined  the  deviation  in  the  angular  spacing  of  the  “spokes”  in  the  distal  shield  of  the  individual 
eoecoliths  provided  in  Fig.  2  of  Gordon  (2006).  For  this  particular  coccolith,  cr^zla  ~  0.27,  and  there 
were  40  open  angular  sectors.  This  coccolith  shield  is  similar  in  size  and  shape  to  the  larger  (2.75  pm) 
pinwheel  (n  =  6,  32  open  sectors)  I  examined.  The  computations  for  the  that  pinwheel  (not  shown) 
suggested  that,  for  the  purpose  of  computing  baekseattering,  the  periodic  pinwheel  is  a  good 
approximation  to  aperiodic  pinwheel  as  long  as  cr^zla  <  0.4  {s  <  1).  This  suggests  that  replacing  the 
aperiodic  fine  structure  of  the  distal  shield  of  E.  huxleyi  eoecoliths  with  a  strictly  periodic  fine  structure 
will  not  degrade  the  modeling  of  their  baekseattering,  especially  for  natural  samples  containing  large 
numbers  of  eoecoliths. 
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